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NOTATION 

Amplitude  of  solitary  wave 

Constants  in  tne  solitary  wave  equation 
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Complex  potential 
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Complex  velocity 

Horizontal  and  vertical  inner  coordinates 
(x  -  x'A',  y  «  y ’ A ' ) 

Outer  coordinates 

Complex  variable 

Dimensionless  wave  number 
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ABS  TRACT 

The  two-dimensional  free-surface  flow  generated  by  a  singu¬ 
larity  moving  with  near-critical  speed  (i.e.  with  Proude  number 
referred  to  the  water  depth  near  unity)is  solved  by  using  the 
method  of  matched  asymptotic  expansions.  In  the  vicinity  of  the 
singularity  the  problem  is  solved  by  an  infinitesimal  wave  ex¬ 
pansion  (inner  emcansion)  while  at  large  distances  from  the  singu¬ 
larity  shallow  water  theory  provides  the  proper  solution  (outer 
expansion).  A  composite  expansion  provides  a  uniformly-valid  so¬ 
lution  for  the  velocity  components,  the  free-surface  profile  and 
the  drag  force.  These  two  basic  approaches  of  water-wave  theory 
are  shown  to  be  solutions  of  the  same  flow  problem,  but  valid  in 
differ  regions.  Although  the  inner  expansion  satisfies  linear 
equations,  the  solution  depends  nonlinearly  on  the  small  parameter 
of  the  problem  (the  singularity  strength). 

INTRODUCTION 

A  singularity  moving  at  constant  speed  in  finite-depth  water 
is  considered  herein.  The  flow  is  two-dimensional  and  it  is 
cteady  when  referred  to  a  moving  coordinate  system.  In  this  sys¬ 
tem  the  standing  singularity  is  considered  to  be  a  perturbation 
of  a  stream  of  uniform  velocity. 

The  nonlinear  free-surface  problem  has  been  solved  by  a 
first  order  infinitesimal  wave  expansion  for  different  types  of 
singularities  (Wehausen  and  Laitone,  i960,  review  the  solutions). 
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All  solutions  diverge,  however,  when  the  unperturbed  velocity 
approaches  its  critical  value,  i.e.  when  the  Froude  number  based 
on  stream  depth  tends  to  unity.  It  is  generally  assumed  that 
the  linear  theory  fails  in  this  case  (Stoker,  1957,  p.  217). 

The  other  basic  approach,  nonlinear  shallow  water  theory,  is 
able  to  predict  the  existence  of  waves  of  finite  amplitude  at 
near-critical  speeds;  but  this  theory  cannot  represent  flows  in 
the  vicinity  of  ^-insularities,  because  it  is  based  on  the  assump¬ 
tion  that  the  horizontal  velocity  component  is  almost  uniform  and 
much  larger  than  the  vertical  one. 

The  purpose  of  this  paper  is  to  solve  the  problem  of  criti¬ 
cal  flow  by  using  the  method  of  matched  asymptotic  expansions 
(Van  Dyke,  1964).  The  infinitesimal  wave  expansion  will  be  shown 
to  provide  a  meaningful  solution  in  the  vicinity  of  the  singularity 
and  will  be  called,  consequently  the  inner  expansion.  This  expan¬ 
sion  diverges  at  a  large  distance  from  the  singularity.  There,  an 
outer  expansion,  derived  according  to  matching  requirements,  de¬ 
scribed  the  flow  conditions.  This  outer  expansion  is  precisely 
the  well-known  shallow  water  theory.  A  composite  expansion, 
therefore,  solves  the  paradoxal  problem  of  critical  flows.  The 
results,  however,  go  beyond  this  purpose;  they  show  that  the  two 
basic  theories  of  water  waves  are  organically  interrelated  by 
the  method  of  matched  asymptotic  expansions. 

The  present  paper  is  an  extension  of  previous  work  in  which 
the  method  of  matched  asymptotic  expansions  was  applied  to  free- 
surface  flows  in  porous  media  (Dagan,  1967). 
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THE  HOMOGENEOUS  PROBLEM:  THE  SOLITARY  WAVE 

The  derivation  of  the  solitary  wave  by  the  method  of  matched 
asymptotic  expansions  will  be  discussed  in  the  first  stage. 

With  the  variables  of  Figure  1  made  dimensionless  by  re¬ 
ferring  the  velocity  components  u,  v  to  (gi1)^  and  the  lengths 
to  t' ,  the  exact  equations  of  flow  are 


u  +  v  =  0 
X  y 


u  -  v  =  0 

y  x 


~\ 


(o  <  y  <  n) 


[1] 

[2] 


i(ua  +  V2)  +  T) 


V 


)  (y  -  n) 


[3] 

[4] 


v  *  0 


(y  «  0) 


l5) 


n  -  i 


(x  -  0) 


[6] 


•1  *  h 


(|x|  -  ») 


[7] 


uA  being  the  velocity  at  the  crest. 

The  solution  for  the  flow  in  the  vicinity  of  crest  A  is 
sought  by  a  first-order  infinitesimal  wave  expansion  whose  zero 
order  term  is  a  uniform  flow  of  velocity  uq  and  depth  t)  ■  1, 
i .  e. 
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u  =  u  +  eui 
o 


v  =  evi 


T)  =  1  +  €T]i 


u.  =  u-  +  eu. 

A  o  Ai 


[8] 


g  being  the  small  parameter. 

By  substituting  the  expansion  [8]  into  Equations  [l]-[6], 
and  discarding  the  condition  at  infinity  [7]>  one  obtains  the 
usual  linearized  equations 


Ulx  +  vly  =  0  1 

/ 

>  (0  <  y  <  1) 

[9] 

u  -  V,  a  0  | 

iy  ix  ) 

[10] 

U  Ui  +  T)X  =  u  UA 

o  *  '  0  Ai 

>  (y  -  i) 

(ill 

vx  -  u  T»  «  0 

o’lx  ^ 

I 

[12] 

vx  «  0 

(y  -  o) 

[131 

t)i  •  0 

(x  -  c) 

[14] 

Considering,  for  the  sake  of  simplicity,  only  the  region 
x  £  0,  the  elementary  solution  of  Equations  (9)-[  1  may  be 
written  as 
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Ui  =  c  {cosh  a  -  cos  ax  cosh  ay)  +  c  sin  ax  cosh  ay  +  u 
S  0  i 


Vx  =  -c  sin  ax  sinh  ay  -  c  slnh  ay  cos  ax 
s  c 


rii  =  -c  u  cosh  a  (l-cos  ax)  -  c  u  cosh  a  sin  ax 
so  '  'co 


where  a  Is  the  root  of  the  equation 


tanh  a  2 

-  _  u 

a  o 


[15] 

[16] 
[17] 


[18] 


and  u  <  1. 
o 

The  solution  is  continued  in  the  region  x  <  0  by  replacing 
x  by  -x  in  Equations  [15] -[173*  The  two  arbitrary  constants  c_ 

b 

and  c  multiply,  therefore,  a  smooth  solution  at  the  origin  and 
c 

a  eusped  one,  respectively  (equation  17). 


The  above  well-known  solution  is  periodic  and  cannot  satisfy 
the  requirements  of  nonperiodic! ty  and  uniformity  at  |x|  =  °°  ex¬ 
pressed  by  Figure  1  and  Equation  [7].  In  the  vicinity  of  the 
origin,  however,  and  for  small  a  (i.e.  uq  close  to  the  critical 
speed  u  =  l)  the  solution  [ 15 3  - [ 18 ]  may  be  expanded  in  a  power 
series  as  follows 


u 


f  € 


c  ax 

c 


+  c. 


a2 (x2-y2+l ) 


+  u 


Ax 


u 


[19] 
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rj  =  1  +  e 


tanh  a 
a 


-  cga2xy) 

[20] 

2  2  \ 

a  x 

c  u  ax  -  c  u  — — — 

[21] 

co  s  0  2  1 

2 

a  .  2 

3  0 

[22] 

with  the  highest  order  retained  being  0(ea2) 


The  solution  expressed  by  Equations  [l9]-[22]  becomes  un¬ 
bounded  as  x  -*  °°,  the  expansion  being  therefore  singular  there. 
Considering  this  as  an  inner  expansion,  an  outer  expansion  is 
sought  according  to  the  usual  procedure  (Van  Dyke,  1964)  by 
adopting  the  following  outer  variables 


I  l 

X  =  e2x  Y  =  y  N  =  T)  U  =  u  V*  e2v 

[231 

The  outer  variables  of  [23]  have  been  selected  so  that  for 
large  x  the  quadratic  term  in  Equation  [21]  becomes  of  the  order 
of  magnitude  0(l)  with  respect  to  e. 

The  inner  variables  may  now  be  expanded  tentatively  in  an  e 
power  series 

U  =  UQ  +  eUi  +  ... 

V  =  eVx  +  ... 

N  =  N  +  eNi  +  . . . 

0 


[24] 

[25] 

[26] 
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The  Equations  [23l-[26]  represent  precisely  the  shallow 
water  expansion  of  Priederichs  and  Keller  (Priederichs,  19^8). 

The  outer  expansion  of  the  exact  Equations  [.!]-[ 6]  has  been 
carried  out  systematically  by  Laitone  (Laitone,  i960)  and  will 
not  be  repeated  here.  The  first  order  equations  are  the  well- 
known  shallow  water  equations,  the  solution  being  (Laitone,  i960) 

i 

U  =  N  2  =  const.  [27] 

00 

Ui  =  f(X)  Ni  =  -[UQf(X)  +  C]  Vi  =  0  [28] 

The  function  f(X)  satisfies  the  differential  equation 


fYY  -  -2—  f2  -  f  =  0  [291 

AA  5  u  £> 

o  0 


where  C  is  an  arbitrary  constant. 


Since  the  solution  has  to  be  used  only  in  the  outer  zone, 
the  nonlinear  Equation  [29]  may  be  easily  solved  by  iterations 
in  the  vicinity  of  X  =  °o  (with  f(«)  =  0);  but  we  may  take  ad¬ 
vantage  of  the  fact  that  the  exact  solution  of  [29]  is  well 
known  being  in  the  nonperiodical  case  (Laitone,  1963) 


f(X)  =  -  ~  sech2 
0 


i 

(x  +  0) 


[30] 
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where  £  is  an  arbitrary  constant. 

The  constants  u  ,  u.  ,  e,  c  ,  c  ,  C,  U  ,  N  and  8  which  ap- 
O  A  i  SC  0  0 

pear  in  the  inner  and  outer  solutions  have  to  be  found  by  matching 
and  by  using  the  Equations  [7]  and  [27]-  The  matching  is  carried 
out  by  requiring  that  the  inner  limit  of  the  outer  representation 
should  be  equal  to  the  outer  limit  of  the  inner  representation 
(Van  Dyke,  1964). 

The  inner  limit  of  the  outer  representation  is  obtained  by 
replacing  the  outer  variables  in  Equations  [27]~[30]  by  the  inner 
variables  of  Equation  [23]  and  expanding  as  e -*•  0  and  X  is  fixed. 
The  result  for  the  function  f  is 


f  ■  - 


U  k‘ 
o 


1  -  2BDe2x  -  D2(l-3B2)  exs 


+  . .  . 


[31] 


where,  for,  briefness,  D  = 


1  JL 

4  6 

.  U  , 

\  o  / 


1 

2 


,  A  =  cosh  D/3  and  B  =  tanh  D/3 . 


In  Equation  [31]  only  the  terms  necessary  for  matching  at  the 
considered  order  have  been  retained. 

The  substitution  of  [31]  in  Equations  [23] -[28],  rewritten 
in  inner  variables,  yields 


N  -  N  -  cCB2  -  m  e3/2x  -  2^1=2 Bil 


[32] 


u  =  u 


.  +  2CBD  €3/2x  +  CD2(1-3B2) 

2 


0  IU!  U  A* 
0  o 


U  A' 
0 


133] 
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The  comparison  of  Equations  [32]  and  [33]  with  Equations  [193 
and  [21]  shows  that  the  matching  at  zero  order  requires 


N 


o 


=  1 


[3^] 


which  means  that  the  unperturbed  flow  must  be  near  critical.  More- 

1 

over,  a  has  to  be  of  order  e2  in  order  to  make  the  matching  pos¬ 
sible  at  higher  order.  For  the  sake  of  simplicity  we  can  take 

a2c  =  e  and  express  the  inner  expansion  [ 19]  and  [21]  as 
s 

t]  —  1  -  c  e^/2x  -  ^-x2  +...  [35] 

Is  C. 

u  =  1  +  eu.  +  c  e^2x  +  ^  e2x2  +. .  .  [36] 

Ai  C  d 

Since  terms  of  order  e2y2  and  €2  of  u  appear  in  e2U2  only 
and  terms  of  order  e2y  of  v  appear  in  e2V2  (Laitone,  i960),  they 
have  been  discarded  in  Equations  [353  and  [36]. 

It  can  be  easily  checked  that  the  pair  of  functions 


which  appear  in  the  inner  expansion  [ 193  and  [20]  are  exact  solu¬ 
tions  of  the  linearised  Equations  [9l-(13]  when  uq  «  1.  They  are 
the  elementary  first  order  solutions  of  the  homogeneous  linearized 
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equations  of  a  free-surface  flow  at  critical  speed,  counterpart 
of  the  sin  and  cos  solutions  [  15 1  and  [  16]  of  a  subcritical  flow. 
The  Equation  [37]  may  be  used  directly  as  first  order  solutions 
of  Equations  [ 9l -[ 13 1  with  uq  =  1.  The  derivation  by  the  limiting 
process  discussed  here  is  useful,  however,  for  the  case  of  a  singu¬ 
larity  as  well  as  for  higher  order  expansions. 

The  matching  of  the  outer  terms  [32]  and  [33]  with  the  inner 
terms  [35]  and  [36]  is  a  matter  of  algebraic  computation. 


In  the  case  of  the  solitary  wave  a  solution  which  is  smooth 

at  the  origin  is  sought.  Hence,  c  =0  and  the  corresponding 

0 

outer  term  2CBD/QoA2  =0,  i.e.  B  =  0,  A  *  1,  and  £  =  0. 

The  inner  and  outer  expansions,  written  in  inner  variables, 
become  in  this  case 


CZ 

1  -  1  -  T 


U  =  U  +  €U.  +  -77-  X2 

o  Ai  2 


*\ 


N  =  N 


U  »  U  - 


3  C  2  2 
0 


l£  ,  2.  C  2  a 
U  4  TT  7  x 
0  U  7 
o 


J 


[38] 
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The  matching  at  first  order  with  the  supplementary  conditions 

U  =  N  2  and  N  =  H  (for  x  -*  «)  yields 
oo  o' 

i 

Nq  =  Uq  =  uq  =  1  eC  =  1-H  C  =  ( 2/3 ) 2  euAi  =  1-H  [39] 

The  complete  solution  of  the  solitary  wave  should  be  ex¬ 
pressed  by  a  composite  expansion,  but  this  is  not  necessary  since 
the  outer  expansion  is  regular  in  the  inner  zone  and  coincides 
with  the  inner  expansion,  i.e.  the  overlapping  zone  coincides  with 
the  inner  zone.  Hence,  the  whole  solution  is  expressed  by  Equa¬ 
tions  [26],  [28],  [30]  and  [39]  in  inner  variables  as 


[40] 


which  is  the  well  known  solitary  wave  solution.  The  velocities 
at  the  crest  and  at  infinite  are 


UA  “  1  -  EUAl  ■  h 

[41] 


or,  in  dimensional  variables 


which  is,  again,  a  classical  result. 


HYDRONAUTICS,  Incorporated 


-12- 

It  seems,  therefore,  that  the  method  of  matched  asymptotic 
expansions  does  not  provide  any  improvement  in  the  accuracy  of 
the  shallow  water  solution.  This  feature  is,  however,  charac¬ 
teristic  for  the  solitary  wave  which  has  a  uniform  structure  in 
the  vicinity  of  the  crest,  but  will  be  different  in  the  case  of 
a  moving  singularity. 

The  derivation  of  the  enoidal  waves  is  entirely  similar  to 
that  of  the  solitary  wave  and  will  not  be  considered  here. 

The  possibility  of  matching  the  outer  expansion  with  an  inner 

solution  which  has  a  non-zero  derivative  at  the  origin  (c  ^  0) 

0 

is  the  key  to  the  solution  of  flow  past  singularities. 

FLOW  PAST  A  SINGULARITY 

A  submerged  vortex  in  a  uniform,  unperturbed  stream  will  be 
considered  for  the  sake  of  definiteness.  Any  other  type  of  singu¬ 
larity  may  be  treated  in  a  similar  way  (Figure  2). 

It  is  worthwhile  for  the  following  derivations  to  discuss 
first  the  radiation  condition.  In  both  subcritical  and  super¬ 
critical  flows  the  flow  is  unperturbed  at  x  »  -*  (Stoker,  1957)* 
Adopting  the  same  condition  in  the  case  of  a  near-critical  flow, 
it  is  convenient  to  take,  R0  *  1  and  uQ  in  the  inner  expansion 
as  the  unperturbed  depth  and  velocity  at  x  »  *»  rather  than  at 
x  «  0.  Since  the  infinitesimal  wave  expansion  will  be  shown  to 
be  regular,  under  the  above  radiation  condition,  in  the  region 
x  <  0,  the  expansion  may  still  be  called  the  "inner"  expansion. 
Hence,  the  first  order  linearized  expansion, similar  to  Equation 
(8)  is 
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u  =  u  +  6(e)ui  'N 
o 

V  =  B(€)Vi  > 

»  1  +  &(Oni 


[42] 


The  gauge  function  6(e)  represents  the  vortex  strength 
(6(e)  -*  0  as  e  -*■  0).  Its  form  will  be  found  from  the  require¬ 
ments  of  matching  with  the  already  found  outer  solution  [32] 
and  [33l-  One  could  obviously  start  the  inner  expansion  with  an 
e  power  series  but  an  unknown  gauge  function  has  to  be  used  in 
this  case  in  the  outer  expansion  in  order  to  make  the  matching 
possible. 

The  variables  of  the  inner  expansion  [42]  satisfy  the  Equa¬ 
tions  [9],  [10],  [12]  and  [13],  Equations  [ 11 ]  and  [14  being  re¬ 
placed  by 


UQUX  +  Hi  -  0 

(y  *  1) 

[43] 

T)l  *  0 

(x  -  -«) 

144] 

As  usual  the  velocity  components  are  separated  in  regular 
and  singular  parts 


[451 


r  .  s 

Ux  «  Ui  +  UX 


r  s 

Vi  -  Vi  +  Vi 
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The  singular  term  is  selected  to  represent  the  potential 
flow  due  to  a  vortex  located  at  x  =  0,  y  =  h  and  confined  between 
two  solid  boundaries  y  *  0  and  y  =  1  (Figure  2). 

The  complex  potential  of  this  flow  is 


which  is  obtained  from  Equations  (43]  and  [12)  by  the  elimination 
of  tjx  •  The  solution  may  be  represented  by  integral  Fourier 
transforms  (Vfehausen  and  Lai  tone,  i960)  as 
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wr  (2)  =  ui 
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[491 


where  u^2  <  1.  Equation  [^9]  will  be  integrated  for  uq  <  1  and 

u  will  tend  to  Its  critical  value  u  =  1  in  the  final  result, 
o  0 

after  Imposing  the  radiation  condition. 

The  details  of  the  integration  for  ui  (x,l)  are  given  In  tne 
Appendix.  The  final  result  for  ui(x,l)  is 


ux(x,l)  *  uiA  ( x ,  1 )  +  ui  (x,  1 ) 
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for  x  2  0.  For  x  f  0,  x  has  to  be  replaced  by  -x  in  Equation  (50], 
c  is  an  arbitrary  constant.  In  order  to  ensure  the  radiation  re¬ 
quirement,  the  wrst  term  is  removed  from  Equation  (50)  for 
x  \  0  by  adding  to  u4  the  proper  solution  of  the  homogeneous  equa¬ 
tion  ( **Q )  (Lamb,  1 9-* 5 ) *  which  Is,  in  the  present  case 
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The  limiting  case  of  u  1  may  be  computed  by  taking  into 
account  that  kQ  is  a  root  of  the  following  equation 


tanh  k  k-k  /3+. . 
k  k 


151] 


Hence  uq2  =  1  -  kQ2/3  +  ...  and  the  final  expressions  of  ui  (x.l) 

when  u  -*•  1  become 
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[53] 


the  higher  terms(in  kQ2,  k*  . . .  )  being  neglected. 

The  complete  solution  of  tj  ,  according  to  Equations  [42], 
(43],  f 1,  (52]  and  [53]  is 
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where  tan  k  =  k  and  k  >  0. 
mm  m 

The  inner  solution  Is  singular  as  x  =  <».  Using  the  previous 
outer  expansion  (Equation  23),  the  inner  solution  may  be  matched 
with  the  outer  solution  [32]  and  [ 33 J  in  the  region  x  >  0. 

The  outer  limit  of  the  inner  representation  is  found  from 
Equation  [55],  through  the  usual  procedure,  to  be 


q  *  1  -  3b6(e)x  [56] 

-kme- 5  X 

where  terms  containing  e  have  been  discarded  since  they 

are  exponentially  small. 

The  outer  limit  of  the  inner  representation  of  u(x,y)  should 
be  determined  by  integrating  Equation  [46]  for  arbitrary  y.  It 
is  known,  however,  that  at  x  >  1  the  u  distribution  is  almost 
uniform  and  u(x,l)  may  be  substituted  for  u(x,y).  Hence,  ac¬ 
cording  to  Equations  [42]  and  [54] 
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[57] 


and  the  outer  limit  is 


u  =  +  3b6x. 

o 


[56] 
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The  matching  of  the  inner  solution  [56]  and  [58]  with  the 
outer  solution  [32]  and  [33]  at  zero  order  implies  N  =  U  =  1. 
The  matching  at  e  order  requires  that  Uq  and  Nq  should  be  e  de¬ 
pendent  (a  similar  situation  occurs  in  the  problem  of  inviscid 
flow  past  an  airfoil  (Van  Dyke,  1964)).  Hence, 


N  =  1  +  eCB2 
o 


U 

0 


u 

0 


[58a] 


These  two  relationships  and  Equation  [27]  show  that,  at  the  order 
considered  here. 


u  =  1  +  eG 


i  B2  -  1 


[59] 


where  the  relationship  1/A2  =  1-B2  has  been  used. 

The  matching  of  the  terms  dependent  on  x  provides  the  ad¬ 
ditional  relationship 


Ce 


[60] 


From  Equations  [591  and  [60]  we  finally  find 
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[61] 
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eC  = 


33(1 -B2  ) 


Equations  [6l]  and  [62]  permit  the  determination  of  B  (also  A 

and  8)  and  eC  as  functions  of  6  and  1-u  . 
r  o 

Two  extreme  cases  are  of  particular  interest:  In  a  first 

case  assume  that  6  -*  0,  but  1-u  ^  0.  Three  solutions  are  then 

o 

possible:  (i)  B  =  1,  i.e.  /3  -  »,  which  means  a  uniform  flow; 

(li)  B  =  0  and  6/B  =  0  which  is  again  a  uniform  flow  since 
eC  =  0;  (iii)  B  =  0  and  6/B  /  0,  which  represents  a  solitary- 
wave.  The  last  case  is  nontrivial  but  has  to  be  excluded  from 
physical  considerations.  It  has  been  considered  as  a  mathemati¬ 
cal  possibility  by  Filippov  (i960)  using  other  methods. 

In  the  second  case  uq  =  1  and-  §  /  0,  which  means  that  the 
unperturbed  upstream  flow  is  exactly  critical.  In  this  case 


1/3  2/3 
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The  outer  solution  for  the  free-surface  profile  in  the 
general  case  according  to  Equations  [26],  [28]  and  [30],  is 


given  by 
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All  the  parameters  appearing  In  Equation  [64]  may  be  determined 
as  1‘unctlons  of  (1-uq)  and  6  from  Equations  [58]*  [ 59 ]  >  [6l]  and 
[62].  A  composite  expansion  with  a  uniform  representation  of  q 
may  be  written  by  combining  the  inner  and  outer  solutions  (Van 
Dyke,  1964).  In  the  simple  case  u  =  1  the  composite  expansion 
of  q,  for  x  >  0,  is 


In  Figure  3  the  shape  of  the  free-surface  is  represented  for  the 
particular  case  b  =  0.6  and  6  =  0.01.  The  series  in  Equation 
[65]  has  been  summed  by  using  an  electronic  computer  with  the 
values  of  from  Carslaw  and  Jaeger  (1959)-  The  matching  is 
possible,  at  the  order  considered  here,  only  for  6  «  1. 

The  drag  force  on  the  singularity  may  be  found  by  a  mo¬ 
mentum  balance  between  the  sections  x  =  -  »  and  x  =  <».  The  flow 
depth  and  velocity  at  x  *  +00  (downstream)  are 
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The  dimensionless  drag  forces  is,  therefore,  given  at  e 


order,  by 


F  =  i  -  u  2-  ill-  —  \  lu  +  [l-Z§ 
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At  the  critical  speed  u  =  1,  F  has  the  simple  expression 


=  oil 


F  =  2 


V3  2/3 


[68] 


or  in  a  dimensional  form 
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DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

Three  basic  lengths  are  present  in  the  problems  discussed 
in  the  foregoing  sections:  the  flow  depth  £',  the  amplitude 
(associated  with  e)  and  the  ratio  u'o2/g.  In  the  infinitesimal 
wave  expansion  the  last  two  lengths  are  referred  to  l '  and  the 
periodic  solution  of  the  homogeneous  equations  in  the  subcritl- 
cal  range  provides  the  well-known  dispersion  relationship  be¬ 
tween  Fr  *  uo'/(g l')^  and  the  wave  Jength  L'  (Equation  18).  For 
small  values  of  a,  where  a  *  2irl'/L',  Equation  [l8]  yields 


3(1-Fr* ) 
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As  Fr -♦  1,  L'  -*■  oo  and  the  homogeneous  problem  has  a  non¬ 
periodical  solution  divergent  at  infinity  (Equations  37)-  The 
shallow  water  solution  is  found  by  relating  x'  to  L1  rather 
than  to  i '  and  by  imposing  a  nonlinear  relationship  between  the 

amplitude  and  the  wavelength  (in  the  case  of  the  solitary  wave 

i 

a  ~  e5).  In  this  sense  the  results  obtained  by  using  the  method 
of  matched  asymptotic  expansions  are  related  to  other  derivations 
of  the  first  order  equations  of  shallow  flow  (see,  for  instance, 
Benjamin  (1967),  p.  561-562). 

Although  the  flow  has  been  found  in  the  vicinity  of  the 
singularity  by  linear  equations,  the  dependence  of  the  solution 
on  the  problem's  small  parameter  (l-h  in  the  case  of  the  solitary 
wave  and  6  is  that  of  the  singularity)  is  nonlinear,  even  at  first 
order  (Equations  [40],  [65]  and  [68]).  In  this  sense  the  solu¬ 
tions  may  be  regarded  as  nonlinear. 

The  solutions  for  the  higher  order  terms  of  the  expansions 
and  their  matching  should  provide  very  valuable  information  and 
insight  into  the  problem,  but  the  computations  become  tedious  as 
the  order  is  increased  in  the  inner  expansions. 

The  method  of  matched  asymptotic  expansions  proves  itself 
to  be  a  powerful  too.  in  solving  nonlinear  problems  of  fluid 
mechanics  in  general,  and  of  water  waves  in  particular. 
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APPENDIX 

THE  INTEGRATION  OP  EQUATION  [49] 


One  of  the  integrals  of  [49]  may  be  evaluated  in  the  fol¬ 
lowing  closed  form  (Gradshteyn  and  Ryzhik,  pp.  505) 


cosk  (z-s) 
cosh  7TS  +  cos  7Tb 


ds 


*00 


slnh  kb  cos  kz 
sin  7rb  sinh  k 


[TO] 


Hence,  Equation  [49]  becomes 
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[71] 


On  the  upper  boundary  z  ■  x  +  i  the  real  part  of  [72]  is 
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[73] 


The  Integration  of  Equation  [73]  follows  Lamb's  procedure. 
For  x  >  0  the  Integration  is  carried  out  in  the  upper  k  complex 
plane.  The  poles  of  the  integrand  have  the  following  locations: 
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two  poles  k  =  ±  k  on  the  real  axes,  roots  of  the  equation 

o 

tanh  k  =  uQ2kj  a  row  of  poles  on  the  imaginary  axes  k  *  ikm 
(m  *=  1,2,...)*  roots  of  the  equation  tan  k  =  u^k  and  a  row  of 
poles  k  =  inir  (n  =  1,2,...). 

The  contribution  of  the  last  series  of  poles  is 
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The  same  sum  is  obtained  from  the  integral 


ikx 

slnh  kb  e 
sinh  k 


dk 


which  may  be  integrated  in  a  closed  form  (Gradshteyn  and  Ryzhik, 
p.  507),  the  result  being 
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The  contribution  of  the  poles  k^  is  given  by  the  sum  of 
residues,  i.e. 
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-k  x 

k  sin  k  b  e  m 
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1  -  1  +  k  2u  2  \  sin  k 
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o  / 

The  contribution  of  the  poles  ±  kQ  is  evaluated  by  in¬ 
tegrating  [73]  along  two  semi-circles,  in  order  to  obtain  the 
principal  value.  The  result  is 


k  sinh  k  b  sin  k  x 
0 _ o _ o 

sinh  k  |l  -  —  +  u  2k  2 

°l  uo8  0  0 

These  partial  results  are  summed  in  Equation  [ 50] 

When  x  <  0  the  results  are  similar,  and  are  obtained  from 
those  above  by  replacing  x  by  -x. 
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FIGURI  2  -  A  VORTEX  IN  A  UNIFORM  STREAM 


0  756fi  g-P 


FIGURE  3  -  A  VORTEX  IN  A  CRITICAL  SPEED  STREAM 

(a)  SCHEMATICAL  REPRESENTATION  OF  THE  FREE  SURFACE  PROFILE 

(b)  DETAILED  REPRESENTATION  IN  THE  MATCHING  REGION 
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